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Abstract. We investigate the strong lensing properties of the Sersic profile as an 
alternative to the NFW profile, focusing on applications to lens modelling of clusters. 
Given an underlying Sersic dark matter profile, we study whether an NFW profile can 
provide an acceptable fit to strong lensing constraints in the form of single or multiple 
measured Einstein radii. We conclude that although an NFW profile that fits the 
^ ' lensing constraints can be found in many cases, the derived parameters may be biased. 

& . In particular, we find that for n ~ 2, which corresponds to massive clusters, the mass 

at r2oo of the best fit NFW is overestimated (by a factor of ~ 2) and the concentration 
is very low (c ~ 2). The differences are important enough to warrant the inclusion of 
Sersic profile for future analysis of strong lensing clusters. 
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\Q • 1. Introduction 
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According to ACDM simulations, dark matter is expected to dominate baryonic matter 
on both galactic and cluster scales. Observations, including X-ray studies, gravitational 
lensing and rotational velocity curve measurements, also suggest that the luminous 
matter is only a fraction of the total matter in the Universe, although it dominates 
the dark matter in the innermost regions of galaxies and clusters. Understanding the 
nature of dark matter is one of the key challenges in modern astrophysics. In particular, 
a knowledge of the spatial distribution of the dark matter is crucial for understanding 
its interplay with the baryonic matter and gaining insight into its nature. 

The dark matter distributions for halos from ACDM simulations are usually 
described by the Navarro, Frenk, & White (NFW) profile [I], but recently [2] found 
a better fit using a Sersic profile for the 3D density distribution. That work was 
expanded in [3] which found that the Sersic distribution also provides a good fit to 
the 2D distribution of dark matter halos and found the deprojected Sersic profile to 
give the best fit to the 3D distribution. It is intriguing that the Sersic law, which 
is often fit to the 2D luminosity profiles of elliptical galaxies @J [5j [6j [7J, should also 
describe the surface density profiles of dark matter halos from simulations. In fact, 
work by e.g. [8j suggests that the Sersic profile provides a better fit in particular to the 
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inner regions of dark matter halos, where the DM density profile affects the kinematics 
of the central galaxies in these halos more strongly. As lensing can probe the surface 
density profiles of galaxies (see e.g., P, [[()]) and clusters (see e.g., [TTJ H2]), especially 
in their central regions, it is of interest to compare the lensing properties of the Sersic 
and the NFW profiles; does the use of an NFW profile for lensing mass reconstructions 
introduce a significant bias given an actual Sersic lens profile, or are the two essentially 
indistinguishable for lensing mass reconstructions? 

The question of the slope for the inner profile in lens systems like galaxies and 
clusters has been addressed by several authors before [131 [H] . Others have predicted 
and/or attempted to reconstruct the slope of lensing systems from the observed positions 
and magnifications of multiply-imaged sources [T2l [15]. In this work, we take a very 
theoretically motivated approach: given either an NFW or Sersic profile for dark mass 
distributions, how would their lensing properties differ? We thereby neglect the fact 
that in actual lens systems, the lensing mass also includes a baryonic component that 
is parameterised in a different way. Even if the stellar mass component of early-type 
galaxies can also be described by a Sersic profile, the scale length and central densities 
are very different than those found in simulations for the Sersic profiles for the dark 
matter component. The work presented here therefore only describes the difference in 
lensing properties between NFW and Sersic profiles if the total mass distribution can 
be described by a single Sersic profile. For cluster systems, the dark matter component 
is a more significant contributor to the lensing properties and the baryonic component 
plays a smaller role than for galaxy lens systems. In this sense the results in our paper 
are more readily compared with current lensing constraints on the dark matter profiles 
of clusters. The analytical lensing relations derived in Section 2 are, however applicable 
to all lens systems that can be described by a Sersic profile, irrespective of the scale and 
whether the matter is baryonic or dark. 

Throughout this paper we assume standard ACDM cosmology with Q\ = 0.70, 
Q m = 0.30 and H = 70 km s -1 Mpc -1 {h = 0.7). For numerical lensing calculations we 
place the lens at a redshift of z = 0.3 (a realistic redshift for both galaxy and cluster 
lenses) and the source at z — 10 (chosen to approximate a source at z = oo, although we 
note that for any z > 1 the results are similar). For the given cosmology, 1" corresponds 
to 4.45 kpc at z = 0.3. The paper is organised as follows: We start by introducing 
the NFW and the Sersic profiles in the remainder of Section 1. In Section 2 we derive 
and discuss the magnification properties of the Sersic profile. In Section 3 we study 
and discuss the differences in the lensing properties of a Sersic surface matter density 
and a corresponding best fitting NFW. We study whether applying an NFW fit to an 
underlying Sersic profile leads to a bias in the mass and concentration parameter and 
discuss the implications our results could have on modelling of observed lenses. We 
summarise our conclusions in Section HI 
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1.1. The NFW profile 

The NFW profile [U [16] has been extensively used to fit ACDM halos from simulations 
and in 3D is given by: 

_ 5 c p c 

Pnfw (r/r s )(l + r/r s ) 2 [) 
where r s = r 2 oo/c2oo is the (3D) scale radius, c 20 o is a dimensionless number called 
the concentration parameter, p c is the critical (3D) density at the redshift of the halo, 
r2oo is the radius inside which the mass density of the halo equals 200p c and S c is the 
characteristic over-density for the halo given by 

_ 200 C 200 /gN 

3 m(l + c 20 o) - e 2 oo/(l + c 20 o) ' 
The 2D projection of the NFW profile is frequently used in the modelling of 
gravitational lenses [TTJ [TBI E21 EES]- It is found by integrating equation (jTJ along the 
line of sight, giving 



1 - ^arctanh^ j (X < 1) 
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and the shear is given by 
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where X = R/r s [20]. These lensing relations of the NFW have previously been studied 
and can also be found in e.g. [2U EQl EH]- 



1.2. Sersic profile 

The Sersic law for surface density profiles is given by 




where E ser is the 2D density, R is the 2D radius, n is the Sersic index, b n is a constant 
chosen such that R e is the radius containing one-half of the projected mass and S e is 
the density at R e . The constant b n is found by solving the equation T(2n, b n ) = r(2n)/2 
and can be approximated as b n w 2n - 1/3 + 4/(405n) + 46/(25515n 2 ) [22]. With 
n = 4 the Sersic profile reduces to the de Vaucouleurs profile, whereas n = 1 gives the 
exponential law. A previous study of the lensing properties of the Sersic profile focussed 
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on n = 3,4,5[23]. In [5] the mass distribution of dwarf- and galaxy-sized halos has a 
mean of n ~ 3.0 but n ~ 2.4 for cluster-sized halos with values as low as n = 2. 



2. Lensing and the Sersic profile 

We begin by calculating the total magnification, fi, for the Sersic profile. From standard 



equations for gravitational lensing [24J we have: 



^ = n \~2 2 ( 6 ) 

(1 — k) z — 7^ 

where k = S/S cri t is the convergence and 

c 2 D s 

Scrit = 4ttG (7) 
is the critical mass density, which depends on the angular diameter distances to the 
source (D s ), the lens (D{) and between the source and the lens (As), c is the speed of 
light in vacuum and G is the gravitational constant. For a spherically symmetric mass 
profile, the shear 7 can be calculated as 

7 = — = K ~ K, (8) 

^-"crit 



t(Y) = A ! Y Y>Y,{Y>)dY> (9) 



where 



r 2 

is the mean surface mass density as a function of radius [2] 
Applying this to the Sersic profile, we find 

User = [(1 ^scr)(l + ^ser 2/€ scr )] , (10) 



where 
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eXp [-fo n (-1 + yV»)] , (11) 



crit 



Kser = = — 26 n nexp[6„J — (12) 

^crit * 

and Y = R/R e . We see from equation ( TTOT) that the Sersic profile has two sets of critical 
curves (i.e. curves in the lens plane where the magnification formally goes to infinity), 
which are defined by 1 — R scr = (tangential critical curve) and 1 -t~ K ser — 2ft ser — 
(radial critical curve). 

These equations define two rings and the radii can be found by solving the equations 
numerically. The larger of the rings, given by 1 — R ser = 0, is the so called Einstein ring. 
Its radius is called the Einstein radius, R e i n , which, for the spherical profiles considered 
here, corresponds to the radius at which the mean surface density equals the critical 
surface density. Figure [1] (solid lines) shows a plot of the magnification as a function 
of radius for different values of n and constant E e /E crit and R e . We can see that the 
two critical curves move further apart for steeper profiles, and correspondingly, that the 
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Figure 1. Plot of the magnification properties of the Sersic profile (solid line) and 
the best fitting NFW profile (dotted line) for different values of the Sersic index n, 
and constant E e = 10 8 M Q kpc -2 and R e = 100 kpc. We see that for the Sersic profile, 
the two critical curves move further apart for higher values of n, corresponding to a 
steeper profile. 
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lowest magnification between them decreases. Our result is not in agreement with those 
of [23], who only found tangential critical curves when studying lensing properties of 
the Sersic profile for n = 3,4, 5. We find that this discrepancy is due to a sign error in 
the calculations of |23j and we note that, provided there is a tangential critical curve at 
R = .Rein > 0, there will always be a radial critical curve for the Sersic profile. This is 
because (1 + R ser - 2K sei )\ R=oc = 1 > and (1 + R ser - 2k sct )\ r=0 = 1 - K ser \ R=0 < and 
as both K ser and R ser are continuous functions, there must exist an R > where a radial 
critical curve will occur. 

3. Comparing the Sersic and the NFW profiles 

In this section we compare the lensing properties of the Sersic and the NFW profiles. We 
use the Sersic 2D mass profile as an input, and find a corresponding NFW which gives 
the best reproduction of the lensing constraints, as described below. Following |25j, 
values of R e ~ 100 kpc with 3 < n < 5 are representative parameter values for galaxy 
sized halos, while R e ~ 1000 kpc with 2 < n < 3.5 are a realistic parameter range for the 
Sersic profile for cluster-sized halos, with lower n being found for more massive clusters 
(which are more likely to act as lenses). In particular, we take 500 < R e < 2500 kpc for 
cluster sized halos as input parameters for our simulations. The third input parameter, 
E e , has a minimum value for strong lensing to occur, i.e. S SOI .| fi= o = S e e 6n > S crit . 
Typical Einstein radii for galaxies are of the order of 1" [26J, while for clusters they are 
of the order of 10" with the largest known Einstein radius to date being 45" [TTJ [T2] . 
We will therefore constrain our input Sersic profiles to have Einstein radii in this range, 
and for each given Einstein radii we will vary n and R e , while S e is calculated from the 
three input parameters R e m,n and R e . 

Our goal is to study whether strong lensing data of an accuracy currently available 
can be well fitted by an NFW profile even if the lensing mass distribution actually follows 
a Sersic profile. In particular, we wish to explore whether the corresponding NFW can 
reproduce the strong lensing data, and if so, what the parameters of the NFW profile are 
and how these relate to typical NFW parameter values obtained for simulated cluster 
halos. Therefore we choose to take a Sersic profile as our assumed input and will fit for 
the corresponding NFW profile given certain constraints. 

The constraints from strong lensing are in general the Einstein radius, relative 
locations of the images and their magnification ratio, which are all sensitive to the 
mass distribution within the disc defined by the outermost image. Of these the Einstein 
radius (deduced from the total separation of the multiple images) is the best constrained 
value, while the precise location of images can be affected by nearby substructure and 
the magnification ratio can be biased by both microlensing and dust extinction in the 
lensing galaxy [271 [28l [29] . The inner critical curve may also give additional constraints 
on the mass profile if central images or merging radial arcs are observed. Under the 
assumption of spherical symmetry, time delay measurements between images provide 
a further constraint on the mass in an annulus, bounded by the radii of each of the 
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Figure 2. The top panels show E for the Sersic profile (solid line) for n=2,5 and 
< R < i? c in- Also shown are the best fitting NFW profile (dotted line) as defined 
in Section 13.11 and an alternative fit (dashed line) which is optimised in a tight area 
around i? G in- The lower panels show the corresponding k (upper set of curves) and k 
(lower set of curves) . Although the detailed shape of the NFW varies with the choice 
of fitting, the resulting analysis presented in the paper is not sensitive to it. The 
value of the other input parameters for the Sersic profile were R e = 1000 kpc and 
i?ein = 50 kpc. 

images, but are hard to measure, in particular in strong cluster lensing where most of 
the known lensed images are background galaxies and not quasars. 

We choose two main approaches to the fitting. The first approach assumes that we 
have multiple systems at different redshifts allowing for constraints on the total mass 
within a number of radii. This approach reflects the common situation that several arc 
systems are used to constrain the lensing mass distribution of clusters |30| . The second 
approach assumes that only a single Einstein ring is detected and used for the lensing 
analysis |31j . In this case we fix the Einstein radii of the two profiles to be the same, 
and fit the inner mass density profile as described below in section 3.1. This choice is 
motivated mainly by the fact that the Einstein radius is easy to constrain accurately 
observationally and that the magnification of multiply imaged arcs and arclets in strong 
lensing clusters constrain the surface mass density on a scale given by the Einstein 
radius. 

3.1. Fitting method 

For the case when only one Einstein ring is used to constrain the system, the fit is 
underconstrained, so we must choose an additional constraint on the NFW profile that 
relates it to the input Sersic profile. Since the magnification and shear for a series 
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of observed extended and highly magnified arcs of the same source is sensitive to the 
surface mass density in the inner regions of the cluster, we require that the NFW profile 
fits the input Sersic projected mass distribution in an inner region bounded by R C i n with 
the Einstein radius itself as a constraint. The fitting itself is done numerically, using 
Levenberg-Marquardt least-squares minimisation. As our primary mode of fitting, we 
take the Sersic surface mass density profile at iV = 10000 points, linearly distributed 
from R = to R = R e i n , and applying a constant relative error in each of the points. 
To check the dependence of the results on our precise choice of fitting, we try several 
other fitting methods consistent with the previously described framework. In particular, 
we redo the analysis using fits where we distribute the points differently, where we fit 
log S instead of S, and where we alter the weight of the points, giving more emphasis to 
different regions. For example, instead of distributing the points over < R < R e i n , we 
constrain the fit to a tight area around the Einstein radius (0.95i? e in < R < 1.05i2 e i n ) or 
to an area around the Einstein radius and the inner critical curve of the Sersic profile. In 
addition we changed the absolute weights of the points to be equal, putting an emphasis 
on the inner points (where we have less constraints from lensing). Although the detailed 
shape of the best fitting NFW profile varies depending on the choice of fitting (see Figure 
[2] for a plot of the E, R, k from two fits for n = 2 and 5), and we note that we do not 
constrain whether the profile is very cuspy or not in the very centre as could be done 
using central images as constraint, the analysis presented in the rest of the paper remains 
qualitatively the same. This leads us to conclude that the precise method of fitting we 
use does not affect our conclusions. 

When we have multiple Einstein radii, we have enough constraints to find the best 
fitting NFW without additional constraints. As each Einstein radii gives an estimate of 
the mass of the system within that radius, we fit the mass profile (rather than the mass 
density profile) of the NFW to the mass profile of the input Sersic profile. We assume 
that we have 9 independent systems at distances of 10-170 kpc for a cluster sized lens 
[30]. We set the uncertainty of the mass estimated at each of the points to be 10%, 
corresponding roughly to the level of current observational accuracy. We find that the 
results from the lensing analysis remain qualitatively the same as for our fitting of a 
single Einstein radius. 

The figures in the paper are based on the single Einstein radius constraint, fitting 
£ in the range of < R < R ein using the 1/N weighting, unless otherwise noted. 

3.2. General properties 

Figure [3] shows i? e i n as function of n for different values of R e . The figure shows that 
i?ein is a strong function of R e and n for low values of n, while for the higher range of 
n the values level off. A similar plot may be obtained keeping R e constant and varying 
E e , which shows i? ein increasing with E e . As explained above, our interest in the strong 
lensing properties of the Sersic profile as compared to the NFW profile motivates us 
to either fix the Einstein radius or fit for the total mass at multiple Einstein radii 
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Figure 3. The Einstein radius, i? G i n , as a function of n, for R e = 70, 150, 700, 1500 kpc 
(denoted by increasing line thickness) and constant S e = 10 8 M Q kpc -2 . The Einstein 
radius is strongly dependent on R e , with larger values of R e resulting in a larger 
Einstein radius, and is also dependent on n, with larger n giving larger i? C i n . The n 
dependence is strongest for low n ~ 2-3 but levels off for higher n. A similar plot can 
be obtained by holding R e constant and varying S e , showing i? C i n increasing with E e . 

when making further comparisons. In the following figures the Einstein radius is set 
to -R ein = 50 kpc, corresponding to a cluster sized halo, but all the results presented 
also apply to systems with different Einstein radii (e.g. galaxies) and to systems with 
measurements of several Einstein radii corresponding to multiple sources. 

Fitting an NFW profile to a Sersic lens will not always give a good fit to the data. In 
fact, given that the number of parameters of the NFW is less than that of the Sersic one 
does not expect a good fit in general. Different indices n of the Sersic profile correspond 
to different slopes in the inner mass profile, whereas the slope of the NFW profile is not 
a free parameter. In Fig.HJwe show the 'goodness of fit', the \ 2 normalised with the 
lowest x 2 m the ensemble. As expected, good fits (low values of \ 2 ) are obtained only 
for a small region around a specific value of n. We note, however, that in most lensing 
analysis an NFW profile is assumed and the best-fit parameters are obtained, even if 
the actual \ 2 is high. In such cases, our results suggest that a fit using a Sersic profile 
should be used instead and will give a much better fit if the underlying profile is in fact 
a Sersic. Figure [5] shows how the best-fit NFW scale length varies with the input Sersic 
parameters. Higher values of n give rise to steeper profiles and hence a lower value for 
the scale length r s of the best-fit NFW profile. 

3.3. Magnification and image configuration 

In Figure [1] we overplot the magnification for the best fit NFW profile for n = 2,3,4,5 
and constant R e and S e (note that for this plot we do not keep i? e in constant). The 
figure shows that the difference in the magnification properties is not very large, except 
for n = 2. For n = 3, the fits are very similar, whereas for n = 4, 5 the main difference 
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Figure 4. The reduced x 2 of the best-fit NFW as a function of the Sersic parameters 
n and i? e , normalised with the lowest x 2 in the ensemble. Some set of parameters give 
significantly worse fits than others, consistent with the fact that the number of free 
parameters for the NFW is less than for the Sersic profile. The Einstein radius is kept 
fixed at R cnl = 50 kpc, as described in the text. 
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Figure 5. Scale lengths, log (r s /kpc), for the best-fit NFW profile as a function of 
Sersic parameters, R e and n. The scale length is strongly dependent on n, with low n 
giving rise to high scale lengths. The Einstein radius is fixed as in the previous figures. 

lies in the location of the inner critical lines. The exact location of the inner critical 
lines depends to some extent also on our choice of fitting, and is therefore not a very- 
robust measurement of differences in the two profiles. In addition, the inner critical line 
is usually observationally less well determined than the Einstein radius. Measurements 
of different Einstein radii, due to sources at different redshifts, are a more common and 
practical additional constraint on the inner mass profile. We found that constraining 
the mass profiles in this way at several different Einstein radii results in very similar 
magnification profiles. 

For a more quantitative analysis of the magnification differences, we define r(/i) as 
the ratio the area on the source plane where magnified images (i.e. images with /i > 1) 
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Figure 6. The ratio of the number of images with amplification equal or greater than 
H to the total number of amplified images, r{fj), for the Sersic (solid line) and the 
NFW (dotted line) profiles for different values of n. The calculation of r(/x) is done in 
the source plane, corresponding to magnification maps of size 4i? C j n times 4i? e j n in the 
image plane, shown in Figure[7l with i? C j n = 50 kpc and R e = 1000 kpc (corresponding 
to a cluster sized lensing halo). The plots show that the magnification properties are 
similar for the two profiles, although the NFW over predicts highly magnified images 
for n ~ 2, while it to a smaller degree under-predicts them for n > 3. 
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with magnification greater or equal to /i occur. Taking R ein to be constant, in Figure E] 
we plot r(fi) calculated over an area in the source plane, corresponding to 4i? ein times 
4i? e i n in the image plane (shown for the Sersic profile in Figure [7J) . Figure [6] shows that 
the difference in the magnification is indeed small, although for n = 2 the NFW profile 
is more likely to produce highly magnified images, whereas for higher n the reverse is 
true to a smaller degree. As the difference in the magnification properties of the Sersic 
and the best fitting NFW profiles is small, our results demonstrate that, given current 
observational accuracy, an NFW profile is adequate to describe the magnification for 
these strong lensing systems. 

In many strong gravitational lens systems the observed images do not lie exactly 
on the Einstein radius, but are displaced in an asymmetric fashion; one image being 
within the Einstein radius and another outside. To illustrate the difference between 
Sersic and NFW models we show the image displacements in Fig. [7] for two different 
choices of the source position. For a fixed source position, the image positions for the 
best-fitting NFW profile is indistinguishable from that of the Sersic profile for high n, 
while for low values of n, there is a slight offset between the image positions, in line with 
the findings from above: NFW profiles produce very similar lensing signatures to Sersic 
profiles, and are distinguishable only for low values of the Sersic index n < 3 unless the 
goodness of fit can be used to separate the profiles. 

3.4- Mass and concentration 

Lensing, both strong and weak, is frequently used to estimate the mass and the NFW 
concentration parameter of clusters, which are then compared to those from simulations 
or X-ray measurements, with lensing usually finding a higher mass and concentration 
[TTJ [321 G31 G2B EI] • Therefore, it is of interest to study whether an underlying Sersic 
profile could cause a bias in the mass and c estimate, if an NFW profile is assumed in the 
strong lensing modelling. As strong lensing is only sensitive at small radii, one should 
keep in mind when making these comparison, that a small deviation in the assumed 
profile can lead to large errors in both the mass and concentration parameter. However, 
although weak lensing should be more accurate for such comparisons as it probes a 
similar distance scale, it is very sensitive to contamination of the background galaxy 
catalogue by cluster members. In fact, [30] find for Abell 1689, that the concentration 
parameter derived from strong lensing is consistent between different data sets and 
methods, while that derived from weak lensing still shows large scatter which are an 
artefact of the method applied. In Fig. El we show the ratio of the projected masses 
for the input Sersic and the best-fit NFW within r 20 o (of the NFW profile). The mass 
ratio depends most strongly on the Sersic index n. Roughly speaking, a higher value 
of n gives rise to lower projected masses for the corresponding NFW fits. Therefore, 
applying an NFW profile to a Sersic lens can lead to a significant error in the mass 
estimation, with the mass at r2oo being overestimated by a factor of ~ 2 for low n. 
Similarly, Figure [9] shows the concentration parameter of the best fit NFW profile as a 
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Figure 7. The image plane magnification map for the Sersic profile for different 
values of n with R cin — 50 kpc and R e = 1000 kpc (corresponding to a cluster sized 
lensing halo). The images are centred on the lens and are of size 4i? C j n times 4i?, C i n . 
The black circles correspond to the tangential and radial critical curves for the Sersic 
profile, and the red circles are the critical curves of the best fit NFW profile, where the 
magnification formally goes to infinity. Also shown are the image geometries for double 
images of the Sersic (diamonds) and the NFW (plus signs) profiles for two different 
source positions (blue and green). The symmetrical configuration of images are not 
distinguishable for any n, whereas for an asymmetric configuration the images are only 
marginally distinguishable for n = 2. 

function of the Sersic index n and scale R e . The figure shows how the concentration 
parameters increase mostly with increasing n from very low concentrations to c ~ 17 
for high n. Extreme values of the concentration parameter (i.e. c < 3 or c > 10), can be 
found and correspond to the input Sersic parameters where the NFW fit is worse and 
where the mass estimate is biased (see Figs. H]and[E]). 

Therefore applying an NFW profile to an underlying Sersic profile, can lead to 
wrong mass estimates. For cluster sized halos with n ~ 2, fitting lensing data with an 
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Figure 8. The ratio of projected Sersic input masses over the best-fit NFW masses 
within r2oo as a function of R e and n of the input Sersic profile. For low n the NFW 
profile overestimates the mass by a factor of around ~ 2, while the opposite holds for 
high n. The Einstein radius is kept fixed at i? C i n = 50kpc, as in the previous figure. 
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Figure 9. Concentration parameter for the best-fit NFW halo as a function of the 
input Sersic profile parameters n and R e . For low n, the concentration parameter 
takes on values of c which are lower than those expected from simulations, while for 
higher n the values of c become higher than expected. Therefore, although the lensing 
properties can be reproduced, the parameters of an NFW fit to a Sersic lens can become 
unrealistic (see also Figs. [3] and [5]). The Einstein radius is kept fixed as in the previous 
figure. 

NFW profile would give values of c which would be lower than those expected when 
compared to CDM simulations. Such low values of c have been found in strong lensing 
analysis, although in general lensing overpredicts c compared to X-ray measurements 
[341 l3~T] . For higher values of n, this trend gets reversed: for n > 3 concentration 
parameters of NFW fits are higher than those seen in simulations. However, such high 
values of n are not expected to describe cluster size halos, but galaxy scaled halos [3]. 

Observationally it is thus likely that NFW lens models of clusters that give rise 
to mass or concentrations which seem inconsistent or only marginally consistent with 
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Figure 10. The difference in the projected mass vs. the mass of the Sersic profile 
at radii R > R C i n for n = 2,3,4,5 (denoted by increasing line thickness) with 
^cin = 50 kpc and R e = 1000 kpc. The NFW fit to the Sersic profile, extrapolated 
to higher radii, underestimates the mass for an underlying Sersic profile for n > 3, 
whereas it overestimates it for low n. 



simulations or X-ray measurements would be better described by a Sersic profile. We 
stress that these results are consistent for all the fitting methods we apply, and whether 
we apply one or several sets of Einstein radii as constraints. 

3.5. Extrapolating to the weak lensing regime 

We have focused on comparing the difference in the strong lensing properties of the 
Sersic and the NFW profiles, making R em the most relevant distance scale. A thorough 
analysis of the difference in the weak lensing properties of the two profiles would require 
optimising the fits in the weak lensing regime and is beyond the scope of this paper. 
However, we can compare the predicted total mass and the weak lensing shear 7 in the 
weak lensing regime by extrapolating the derived functions to higher radii. We have 
already compared the projected mass within r2oo for different values of n and R e (see 
Figure [8]), but here we study the NFW: Sersic mass ratio as a function of radius for a 
halo with R e = 1000 kpc and _R e m = 50 kpc, corresponding to a cluster sized halo. In 
Figure [T0l we plot the difference in the projected mass over the projected mass of the 
Sersic profile. We see that for n — 2, the best fitting NFW profile overestimates the 
mass by a factor of > 2 for R > 10i? e i n , while it underestimates it by ~ 2 for n > 3. This 
is due to the fact that in order to fit the low n = 2 Sersic profiles, the estimated r s of the 
NFW (the scale radius at which the density slope increases) is overestimated, making 
the NFW profile fall off less steeply than expected in the outer regions. Similarly, in 
Figure [HI we plot the difference in the projected shear 7 as a function of radius. We 
find that for n > 3, the difference is of the order of 0.01 while for n = 2 it is of the order 
of 0.05. As weak lensing measurements of clusters are usually looking at shear signals 
around 0.01 — 0.05, this is a significant difference, in particular for the low n. Therefore, 
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Figure 11. The difference in the shear of the Sersic and NFW profiles at radii 
R > i?cin for n = 2,3,4,5 (denoted by increasing line thickness) with i? e j n = 50 kpc 
and R e = 1000 kpc. For n > 3 the NFW fit to the Sersic profile, extrapolated to higher 
radii, underestimates the shear by ~ 0.01. For low n, the shear can be more strongly 
affected: it is underestimated by up to ~ 0.05 for low radii and overestimated by a 
similar amount for R > 10i? c in- 

using an NFW to model an underlying Sersic profile with low n would significantly 
over predict the weak lensing signal at large radii, and, depending on the measurement 
errors, result in a poor fit. 

3.6. Observational implications 

We have taken a Sersic profile as the underlying description of the matter profile and 
studied how lensing analysis applying an NFW fit would differ and lead to biased results, 
depending on the input Sersic index. In particular we have found that for n ~ 2 the 
mass can be overestimated by a factor of two and give low concentration. We have 
however not addressed which profile is a better description of real dark matter halos. 
To do so, one should model a real system using the Sersic profile and compare it to the 
results found when using an NFW profile. While such a study is beyond the scope of 
this paper, we discuss here the implications our results have on lens modelling of real 
systems. 

For galaxy scale lenses, various authors [35j EEJ [15] have established that lenses on 
average are well fitted by an isothermal profile, although there can be significant scatter 
for individual lenses [15]. These all apply however to the total matter distribution, both 
luminous and dark, while we have been studying the differences for dark matter only 
profiles. More recent studies have attempted separating the dark and luminous matter 
profiles in lensing of galaxies [374 EE], but in general the constraints from lensing alone 
are not sufficient for separating the two. 

Group or cluster scale lenses, with which we are primarily concerned in this work, 
are more dark matter dominated. This means they are also better suited for an 
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observational test of whether real clusters are better described by an NFW or a Sersic 
profile. Regarding observational strong lensing constraints the main observable is the 
total mass and the slope of the mass profile in the innermost (strong lensing) region of 
the cluster or group. The NFW and the Sersic profile differ in that the Sersic profile 
has an extra independent parameter, n that directly relates to the inner slope of the 
mass distribution. Even though the concentration parameter of the NFW profile is 
related to the slope of the profile, it is also related to the overall scale and mass of the 
cluster and is not an independent parameter. From an observational perspective the 
key quantities then are the slope and mass of the dark matter in the inner region of real 
clusters, and whether they are consistent with an NFW profile as expected from CDM 
simulations. As pointed out recently |3T] NFW fits to observed strong lensing clusters 
very often give masses or concentration parameters that seem inconsistent with CDM 
predictions. In general the trend is that observed clusters have higher concentration 
parameters than predicted in simulations - but some lensing clusters have also yielded 
very low concentration parameters when fit with an NFW profile (e.g. MS1137. 5+6625 
[39] and MS 2053.7-0449 [1Q]). Clusters with unusually high or low concentration values 
may very well be fit much better with a Sersic profile - since the Einstein radius then 
becomes a function of n, even for a fixed scale, R e , an index n (and hence slope) and 
a total mass M may be found that match the predicted shape of the mass profile from 
CDM more closely. It is important to note that the lensing fits of clusters are often 
quite poor, i.e. the \ 2 is high or the rms of the predicted image position is larger than 
the uncertainty in their measured positions [301 14"0] . 

4. Discussion and Conclusions 

In this paper, we have compared the Sersic and NFW dark matter surface density profiles 
with respect to their strong lensing properties. Taking Sersic profiles, with parameters 
that are in the range found for clusters and galaxies in large N-body simulations, we 
explored the parameters of the NFW profile that fits the inner projected surface mass 
density of the Sersic best. The NFW profiles were constrained to have the same Einstein 
radius as the input Sersic profile, as this is the most accurately determined observable 
in lens systems. As an alternative approach, we have assumed a data set consisting 
of several multiply imaged systems at different redshifts, giving a constraint on the 
mass profile at different radii. We find that an NFW profile can in general accurately 
produce the magnification and image positions of a Sersic matter distribution only for 
n > 3. For lower Sersic index values, which are more likely to occur in cluster sized 
halos, the difference in the lensing properties increases, to the level where it could affect 
results of lens modelling, in particular for the magnification estimates. We also note that 
although one can in general find an NFW profile which gives similar lensing properties, 
the parameters of the profile may become unrealistic. For n ~ 2, as simulations suggest 
can be the case for clusters, the mass is overestimated by a factor of ~ 2 while the 
concentration parameter is ~ 1-3, which is lower than expected from simulations, 
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although such low values have been seen in lensing studies [21] • Extrapolating into 
the weak lensing regime we find that the NFW overestimates both the total mass and 
the shear for n = 2, while it underestimates them for n > 3. 

Therefore, if an underlying Sersic profile with n ~ 2, as is found for massive systems, 
is fit with an NFW using strong lensing constraints, it will overpredict both the mass 
(by a factor of a few) and the shear (up to ~ 4%) of the halo at large radii, and 
could therefore contribute to explaining why lensing mass estimates are greater than 
those found by X-rays. Conversely, weak lensing data can in principle be used to 
demonstrate whether the Sersic or the NFW profile is a more accurate description of 
the mass distribution. However, in practice a good S/N ratio is required - even for 
Sersic profiles with low n an accuracy of better than ~ 4% is needed to distinguish the 
profiles and identify the shape of the mass profile accurately enough to obtain reliable 
parameters. Most weak lensing data is not yet of sufficient accuracy to rule out an NFW 
profile as a good representation of the mass distribution even if the mass distribution is 
in fact a Sersic profile with low n. 

In our analysis we have neglected the contribution of baryons - our results are 
therefore most applicable to massive systems, like clusters of galaxies where the dark 
matter dominates. Even though our results clearly show that for a range of Sersic 
parameters, in particular for halos with low n (corresponding to the most massive halos), 
strong lensing data can not always be well reproduced by an NFW profile, we do not 
address the question which of the two mass distribution is in fact a better fit to the mass 
distribution in clusters. However, the results from N-body simulations suggest [3] that 
the Sersic profile is a better fit, at least in the inner regions, which are also the regions 
probed by strong lensing. Even if the overall mass profile of clusters is often fit well by an 
NFW profile, the NFW profile may not reproduce the mass distribution in the inner part 
adequately. As demonstrated here, the sensitivity of the strong lensing properties on 
the form of the mass profile in that region then means that fitting an NFW profile to the 
lensing data will lead to rather meaningless parameters; the concentration parameters 
and masses of NFW profiles recovered in this way can not simply be compared to the 
values obtained in simulations. In this work we have for the first time discussed the 
strong lensing properties of the Sersic profile, given important lensing relations, and 
compared it to the NFW profile. The differences are important enough to warrant the 
inclusion of Sersic profile for future analysis of strong lensing clusters. 
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